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ABSTRACT 

Context. RX J1713.7-3946 is the brightest shell-type Supernova remnant (SNR) of the TeV y-ray sky. Earlier Fermi-LAT results on 
low-energy y-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent 
with a hadronic origin. 

Aims. We update the GeV-band spectra using improved estimates for the diffuse galactic y-ray emission and more than doubled data 
volume. We further investigate the viability of hadronic emission models for RX J1713.7-3946. 

Methods. We produced a high-resolution map of the diffuse Galactic y-ray background corrected for the HI self-absorption and used 
it in the analysis of more than 5 years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport 
equation to calculate the acceleration and propagation of cosmic-rays in SNR. We then determined spectra of hadronic y-ray emission 
from RX J1713.7-3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux 
variations that would allow to test the model with observations. 

Results. We find that RX J1713.7-3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphol¬ 
ogy is best described by that seen in the TeV band. The measured spectrum of RX J1713.7-3946 is hard with index y = 1.53 + 0.07, 
and the integral flux above 500 MeV is F = (5.5 + 1.1) x 10"^ photons cm“^ s"k We demonstrate that scenarios based on hadronic 
emission from the cosmic-ray precursor region are acceptable for RX J1713.7-3946, and we predict a secular flux increase at a few 
hundred GeV at the level of around 15% over 10 years, which may be detectable with the upcoming CTA observatory. 
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1. Introduction 

RX J1713.7-3946, also known as G347.3-0.5, is a young 
(^SNR - 1600 yr) shell-type SNR located in the Galactic 
plane within the tail of the constellation Scorpius. This SNR 
was discovered in 1996 during the ROSAT X-rays all-sky sur¬ 
vey ( Pfeffermann & Aschenbach|1996 ). Its shape is slightly el¬ 
liptical with a maximum extent of 70', and it exhibits bright 
X-ray emission predominantly from the western edge of the 
shell. The assumed distance Jsnr - 1 kpc implies a shell ra¬ 
dius rsNR - 8.7 pc. The remnant contains an X-ray point-like 
source whose properties are similar to central compact objects 
in other SNRs, thus suggesting RX J1713.7-3946 is the remnant 
of a core-collapse supernova. 

Observations with the Japanese ASCA satellite ( [Koyama| 
|et al.||199T I showed a featureless X-ray spectrum in the north¬ 
west shell of the remnant, clearly indicating non-thermal emis¬ 
sion. Two years later, new ASCA observations ( |Slane et al.| 
1999| ), confirmed the absence of line emission everywhere in the 
remnant and showed that all parts of the remnant had power-law 
with indices between ~ 2.2 and 2.4. More recent studies of the 
remnant conduc ted with Chandra ([Uchiyama et al. |2Q03) ) and 
XMM-Newton ( [Cassam-Chenai et al.| 2004| ) satellites clearly 
showed small-scale spatial variations of the photon index rang¬ 
ing from 1.8 to 2.6 which are tracers of the acceleration and 
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propagation history of recently accelerated electrons (Rettig & 
|Pohl|2012| ). 

At higher energies two competing radiation processes com¬ 
plicate the interpretation of the emission from RX J1713.7- 
3946. In the leptonic scenario the y-rays are produced by elec¬ 
trons via inverse Compton scattering, and in the hadronic sce¬ 
nario the y-rays are due to -decay from proton-proton interac¬ 
tions. Both scenarios can produce similar fluxes in the GeV-TeV 
energy range. 

The H.E.S. S. telescope array detec ted RX J1713.7-3946 at 
the TeV-scale ( [Aharonian et al. [2QQ4| ) and showed a similarity 
between the X-ray and TeV-band morphology. No decisive con¬ 
clusion was drawn for the particle population responsible for the 
emission, although the hadronic scenario was favoured. 

The Large Area Telescope (LAT) ( [Atwood et al.||2009| ), the 
principal instrument on board NASA’s Fermi y-ray satellite, 
measured an extended GeV y-ray emission coincident with the 
position of RX J1713.7-3946 ( jAbdo et "aLlpOllj ). The spec¬ 
tral analysis found a very hard spectral photon index (y ^ 
1.5), which was well in agreement with a leptonic scenario. 
Nonetheless, this conclusion did not exclude that protons are ac¬ 
celerated in this SNR. 

In this paper we reanalyze Fermi-LAT data that have 
more than doubled in exposure since the original publica¬ 
tion. Moreover, we develop dedicated tools to produce high- 
resolution background maps that take into account HI self¬ 
absorption. Our background modeling permits determining the 
key parameters of y-ray emission from RX J1713.7-3946, such 
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as spectral index, flux and morphology, with reduced system¬ 
atic uncertainty. Next, we present a model for hadronic y-ray 
emission from the cosmic-ray precursor of the forward shock of 
RX J1713.7-3946, that is in principle viable, and, what is more 
important, we propose a way to verify the model. Based firstly 
on analytic estimates and then on more detailed calculations 
we predict a y-ray flux increase at the level of 15% on a time 
scale of 10 years. We suggest this flux variation be detectable by 
Cherenkov Telescope Array (CTA |Acharya et al.|2013| ). 


2. Analysis and results 

2.1. LAT data 

The LAT is a pair-conversion telescope designed to detect pho¬ 
tons from ~ 20 MeV to 30 0 GeV. A deta i led de scription of the 
LAT instrument is given in [Atwood et al.| ( [2b09| ). The data used 
in this study come from observations over a period of five years, 
from August 8, 2008 to August 13, 2013. Data analysis is per¬ 
formed with public analysis tools developed by the LAT team, 
using the post-launch P7v6 data selection with the appropriate 
instrument response functions. To minimise the spill-over of at¬ 
mospheric y-rays from the Earth limb, a zenith angle cut of 100° 
is applied together with a method of correcting the exposure for 
the zenith cut itself. The method chosen consists in excluding 
time intervals where any part of the region of interest (ROI) is 
beyond the zenith-angle limit. Furthermore, data are not taken 
in consideration while the observatory is transiting the South 
Atlantic Anomaly (SAA) or when the rocking angle exceeds 52°. 

In the analysis only photons with reconstructed energy 
greater than 500 MeV, for which the 68%-containment radius of 
the PSF is narrower than ~ 1.5°, are used. The choice of select¬ 
ing events above 500 MeV is motivated by the earlier published 
spectrum and by the broad PSF at low energy. The broadening 
may in fact affect the analysis by leading to systematic prob¬ 
lems of source confusion in any densely populated region of the 
Galactic plane. 

The data analysis is done using the publicly available Fermi- 
LAT ScienceTools version 9.31.1. In particular, spatial and 
spectral analyses are performed with dedicated scripts based on 
the Python likelihood tools that expand upon the command line 
tools provided with the Fermi Science Tools package. 

Since the analysis is performed over a large number of ob¬ 
served events, the binned likelihood analysis is preferred. The 
region of interest (ROI) is a square region measuring 28° on 
a side and centred on a = 258.39° and d = -39.76° (J2000), 
the nominal position of RX J1713.7-3946. The angular binning 
is chosen to be 0.125° per pixel in stereographic projection to 
match the resolution of the corrected diffuse y-ray background. 
To preserve accuracy in the likelihood analysis 30 energy bins 
are chosen, which allows to accommodate rapid variations in the 
effective area with decreasing energy below ~ 1 GeV 

The y-ray background is modelled in an acceptance region, 
called the source region, that is 5° larger than the ROI. This 
accounts for sources that lie just outside the data region, but 
whose photons spilled into the data set. The background model 
includes 174 point-like sources and 2 extended s ources listed 
in the 2^^ Fermi-LAT catalog ( Nolan et al.||2012 ). The source 
2FGL J1712.4-3941 is not considered in the model because it is 
spatially coincident with RX J1713.7-3946. Furthermore, we in¬ 
clude an isotropic templat^that accounts for extragalactic emis¬ 
sion and residual cosmic-ray contamination. Finally we use a 
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Fig. 1. Counts map of the region of interest (ROI) used in this 
work. Sources of the 2“^ Fermi-LAT catalog are marked by light 
blue crosses. The green dashed line denotes the Galactic plane 
and the red cross represents the nominal position of RX J1713.7- 
3946. During the fit procedure the parameters of the sources in¬ 
side the green circular region of 5° radius are free to vary. 


modified template for galactic diffuse y-ray emission that shall 
be described in the next section. 

Amongst the point-like and extended sources, 105 lie inside 
the ROI and most of them are modelled with power-law energy 
spectra. The spatial and spectral parameters of sources located 
outside a circular region of 5° radius centred at the nominal po¬ 
sition of RX J1713.7-3946 are kept fixed at the values given 
in the catalog. Inside this region the spectral parameters of 7 
sources along with the normalisation of the diffuse components 
are permitted to freely vary. Thus, the final set of free parame¬ 
ters is reduced to 20, a number which allows the fit to converge. 
Fig. shows the ROI with 105 sources of the 2“^ Fermi-LAT 
catalog and the circular region of 5°. 

2.2. The HISA contribution to the diffuse galactic y-ray 
emission 

The likelihood analysis of data from Fermi-LNI requires tem¬ 
plates for an isotropic component, intended to cover extragalac¬ 
tic emission as well as residual instrumental background, and 
diffuse Galactic y-rays. The latter is based on infrared tracers 
and spectral-line information on atomic hydrogen and carbon 
monoxide as tracer of molecular gas. 

While the absolute normalization of the templates is ad¬ 
justable in the likelihood optimization, its angular structure is 
not. If an extended source is to be analyzed, a poor representa¬ 
tion of the diffuse background will lead to systematic errors in 
the reconstructed morphology and possibly a significant offset in 
the derived spectrum of the source. 

A major source of uncertainty is Hi absorption, to which a 
first-order opacity correction is often applied by assuming a uni¬ 
form spin temperature, T^. In terms of the brightness tempera¬ 
ture, Tb, and for low continuum intensity, a correction factor can 
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be derived as 





( 1 ) 


where one typically uses values of in the range 120-150 K. 
While this method has some merit, high-resolution Hi sur¬ 
veys have recently indicated the existence of large variations in 
spin temperature from one cloud to the next. Small-scale self¬ 
absorption features (HIS A, as for H i self-absorption) were iden¬ 
tified that are caused by relatively cold (~ 50 K) gas clouds that 
are located in front of much warmer gas and sources of contin¬ 
uum emission. Their effect on spectral-line data goes beyond a 
simple correction a la Eq.[2on account of both its absorption ef¬ 
fect on the baseline continuum spectrum and the substantial gas 
mass that can be carried by the compact clouds. 

We h ave used data from the Southern Galactic Plane Survey 
(SGPS) ( [McClure-Griffiths et al.||^05] ) to refine the standard 
model of galactic diffuse emissioi^ to a higher intrinsic angu¬ 
lar resolution than is provided in single-dish H i s urveys and to 
correct for HIS A using the algorithm described in [Gibson et al. 
( [2005) ). 


The idea behind Gibson’s method is to iteratively remove 
large-scale spectral and spatial emission structures from H i data 
and to flag the small-scale negative residuals as self-absorption 
features. The additional optical depth imposed by the cold H i ab¬ 
sorber is estimated using the brightness-temperature difference 
between a line of sight through the absorber and other lines of 
sight in the vicinity, assuming most of the continuum signal is 
produced behind the structure. The density and mass of the ab¬ 
sorber are calculated assuming pressure equilibrium. Further de¬ 
tails on the analysis can be found in Gibson’s paper. 

The additional H i column density derived for the absorber is 
then combined with that associated with the absorption correc¬ 
tion to the original line signal. For both we calculate t he loc ation 
on the line of sight using the method of Pohl et al.| ( |200^ . We 
thenjuse the gamma-ray emissivity derived with th e GAFPRQP 
version 54, with best-fit parameters given in Strong et ^ 
to compute the expected y-ray intensity from the addi- 
column density, that we found in the HISA search. 


codq^ ^ 
( | 2010 | , 


pOTOl 
tional H 

Finally, we correct the standard templates of diffuse galactic y- 
ray emission for the additional y-ray signal. 

As high-resolution gas maps are available only for a narrow 
strip along the galactic plane, they are combined with the low- 
resolution maps of the Feiden-Argentina-Bonn survey ( [Kalberla 


et al.|2005) ), and the composite map is sampled with a pixel size 
of 0.125°. At the location of RX J1713.7-3946 we thus have 
an accurate sampling of structure in the diffuse galactic y-ray 
emission on scales commensurate with the angular resolution 
of FAT, whereas off the plane the intrinsic resolution is that of 
the standard galactic-emission template available at the Fermi 
Science Support Center. 

To study the impact of the correction for the HISA on the 
diffuse Galactic emission (DGE) two binned likelihood analyses 
are performed, one with the standard Fermi-FAT DGE template 
and one with the corrected DGE template. The detailed setup 
of the likelihood fit, used for the standard Science Tool gtlike, is 
identical for both analyses. The two likelihood fits are performed 
by assuming as a spatial template for RX J1713.7-3946 a uni¬ 
form disk of 0.5° radius. The fits are performed in the energy 
range spanning from 500 MeV to 300 GeV and indicate that the 
revised model of galactic y-ray emission in the ROI is favoured 
with a significance of ~ 4 cr. 


^ Publicly available at the Fermi Science Support Center 
^ http://galprop.stanford.edu/ 
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Fig. 2. Map of the difference in fit residuals of the revised and 
the standard template of galactic emission. 


Figure shows the difference in the best-fit model maps of 
the two cases, [mi -m 2 ]/ where mi and m 2 refer to the stan¬ 
dard and the revised templates, respectively. To be noted from 
the figure are the localized negative features, that correspond to 
HISA corrections. At the location of RX J1713.7-3946, no sig¬ 
nificant absorption feature was found, and, as we have tested, us¬ 
ing the revised background model changes the best-fit flux from 
RX J1713.7-3946 by only about 1%. 


2.3. Position and spatial extension of RX J1713.7-3946 


After the first 2 years of science operation of Fermi, FAT has es¬ 
tablished a firm detection of y-ray emission from RX J1713.7- 
3946 with statistical significance of ~ 9cr. Now, 3 years later, the 
accumulated exposure has more than doubled, and it is expected 
that systematic uncertainties arising from, e.g., the spatial distri¬ 
bution of emission coincident with the SNR have an impact on 
the reconstructed y-ray spectrum. Before determining the spec¬ 
trum, we therefore need to first find the localization and spatial 
distribution of emission from RX J1713.7-3946. 

To find the center of gravity of events from RX J1713.7- 
3946, we construct a map of test statistic (TS) ( Mattox et al. 
1996), defined as logarithmic improvement in the likelihood 
function L, 


TS = 2 (InL-lnLo) = 2 ln|^j , (2) 

that we achieve by moving a putative point-like source through 
a grid of locations on a square region. The region measures ~ 3° 
on a side and is centered at the nominal position of the SNR. At 
each point of the grid the parameters of the point-like source 
are free to vary and those of all background components are 
fixed at the values found in the previous likelihood analysis (see 
Section[2^. The resolution of the grid is 0.1° which implies 961 
likelihood analyses. The resulting TS map is shown in Fig. 
The envelope of signal with TS > 36 suggests spatially extended 
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Fig. 3. Map of the test statistic (TS) for a point-like source in the 
region around RX J1713.7-3946. The black cross denotes the 
best-fit centroid of a disk model and the yellow circle indicates 
the energy-averaged 68% confinement region of the point-spread 
function. Red circles indicate the position of the second Fermi- 
LAT catalog sources in the background model. Shown in white 
are contours of gamma-ray excess counts based on H.E.S.S. ob¬ 
servations, the levels are 25, 50, and 75. 


emission, because a true point source would appear with local¬ 
ization uncertainty corresponding to the energy-averaged width 
of the point-spread function (PSF, indicated in Fig.j^by the yel¬ 
low circle) divided by the detection significance, i.e. with < 0.2° 
in diameter. Particularly interesting is the local enhancement of 
the signal in the northwest region of the shell with TS > 49, 
which matches well the excess-event density seen with H.E.S.S., 
here indicated with white intensity contours with levels 25, 50, 
and 75. 

To find the best position of the source and hence minimize 
the systematic uncertainties, we then modeled RX J1713.7-3946 
as a uniform disk with a radius of 0.5° instead of a point-like 
source. The best-fit centroid for the disk model is indicated in 
Eigj^and located aia = 258.32° and S = -39.71° in J2000 with 
an error radius of 0.02° at the 68% confidence level, which is 
about 0.07° off the nominal coordinates of the source in direction 
of the Galactic plane. All the following analyses are performed 
using spatial templates centered at this best-fit position. 

To investigate the spatial morphology of the emission as¬ 
sociated with RX J1713.7-3946 a number of spatial templates 
are tested. Eive templates are uniform disks with radii ranging 
from 0.5° to 0.7° in steps of 0.05°. An additional template is 
designed to reflect the TeV-band intensity distribution observed 
with the H.E.S.S. telescope. In Eig.|^we show all these spatial 
templates after being convolved with the Eermi-LAT PSE. Due 
to the broadening of the PSE at low energies the detailed shape 
and size of the six templates can only be distinguished at high 
energies. Table [^summarizes the best-fit parameters in the like¬ 
lihood analysis of each template, performed in the energy range 
from 500 MeV to 300 GeV by optimizing a spectral model of 


H.E.S.S. r=0.50° r=0.55° r=0.60° r=0.65° r=0.70° 



Fig. 4. Spatial templates used to model the emission associated 
with RX J1713.7-3946. The templates are convolved with the 
Eermi-LAT PSE and are shown for three different energies. 


the form 



( 3 ) 


where Nq is the prefactor, Eq the energy scale, and y is the spec¬ 
tral index. 

All models yield a high test statistic (150 < TS < 163), 
about 3 times that achieved in a point-source fit, implying that 
the emission is well resolved by the EAT as an extended re¬ 
gion. Eor uniform-disk templates TS slowly decreases beyond 
a radius of 0.55°, indicating that this value marks the extent of 
GeV-scale emission from RX J1713.7-3946. The H.E.S.S. tem¬ 
plate fits best, and the TS value for the uniform disk of radius 
0.55° is also acceptable (ATS = 2). As the H.E.S.S. template 
does not fit significantly better than the disk of radius 0.5° that 
we used to find the centroid position, we see no reason to devi¬ 
ate from the sequence of analysis steps defined a priori, namely 
first finding the centroid position and then determining the best¬ 
fitting template. We thus use the H.E.S.S. template, centered at 
a = 258.32° and S = -39.71° in J2000, in the following spectral 
analysis. In the end, the vagaries in choosing the spatial template 
and centroid position contribute to the systematic-uncertainty 
margin that we shall discuss below. 


Table 1. Morphological analysis of the y-ray emission associ¬ 
ated with RX J1713.7-3946. 


Morphology 

Flux^ 

Photon index 

TS 

Disk r - 0.50° 

4.982 + 0.916 

1.508 + 0.075 

160 

Disk r = 0.55° 

5.122 + 0.941 

1.509 + 0.075 

161 

Disk r = 0.60° 

5.449 + 1.069 

1.513 + 0.075 

159 

Disk r = 0.65° 

5.547 + 1.052 

1.514 + 0.075 

156 

Disk r = 0.70° 

5.908 + 1.216 

1.521 +0.076 

150 

H.E.S.S. 

5.522+ 1.075 

1.528 + 0.074 

163 


The integral flux from RX J1713.7-3946 is calculated over the 
energy range 500 MeV - 300 GeV and it is given in units of 10"^ 
photons cm"^ s"k 


2.4. Energy spectrum of RX J1713.7-3946 

As first step toward measuring the spectrum of y-ray emission 
from the SNR, we perform a global likelihood analysis using the 
H.E.S.S. template as spatial model for RX J1713.7-3946. The fit 
yields a y = 1.53 + 0.07 and the measured integral photon flux 
above 500 MeV is F = (5.52 + 1.07) x 10“^ photons cm“^ s“^ 
In order to obtain a spectral energy distribution (SED) for the 
SNR, the entire energy range is divided into 9 logarithmically 
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Fig. 5. Broadband y-ray spectrum of RX J1713.7-3946 as mea¬ 
sured here (black points) and by H.E.S.S. (blue points). The two 
bow-tie bands denote systematic (dark green) and total (light 
green) uncertainties of the single-power-law fit, the latter assum¬ 
ing that systematic and statistical uncertainties can be added in 
quadrature. The histogram indicates the TS value for each en¬ 
ergy bin. The dotted line is the threshold for setting a 95% con¬ 
fident level upper limit. 


spaced energy bins. For each individual bin a likelihood fit is 
performed using the H.E.S.S. template as spatial model for RX 
J1713.7-3946 and the spectral model and parameters obtained in 
the global fit (and shown in Table[^. For those energy bins where 
TS< 7 a 95% confidence-level fiux upper limit is calculated. 

To estimate the amplitude of systematic uncertainties, the 
data are re-fitted for several variations of the default model. First, 
the corrected diffuse-background model is replaced by the stan¬ 
dard background model, yielding a 1% variation in the derived 
fiux. Second, the H.E.S.S. template is substituted with reason¬ 
ably fitting uniform-disk templates (radii 0.5°, 0.55°, 0.6°), in¬ 
dicating an uncertainty of ~6% at 20 GeV increasing to 10% 
at 500 MeV and 8% at 300 GeV. For the worst-fitting tem¬ 
plate, the disk of 0.7° radius, we find an offset of 10% at 20 
GeV Furthermore, uncertainties due to instrumental effects are 
also taken into account. The fiux measurement in fact depends 
on the knowledge of the effective collecting area of FAT as a 
function of energy. The systematic error in the effective area 
for the P7SOURCE_V6 event class may be quoted as 10% at 
100 MeV, 5% at 560 MeV, and again 10% at 10 GeV and 
above ( [Ackermann et al.||2012| ). All these contributions com¬ 
bined indicate overall systematic uncertainties at the level of 
about 11% at 500 MeV, 12% at 20 GeV, and 13% at 300 GeV 

The spectral data points are given in Table and the 
resulting energy spectrum of RX J1713.7-3946 is shown in 
Fig. [5] together with previou sly published H.E.S.S. measure¬ 
ments ( [Aharonian et al.|2()TT] ). The two green bands indicate the 
systematic and the total uncertainties in a fit of a single power 
law. The 9 Fermi-FAT data points (black crosses) only include 
statistical errors. In the background of the figure the gray his¬ 
togram represents the distribution of TS values obtained for each 
bin of the Fermi-FAT SED. The dotted line at the bottom of 
the histogram denotes the threshold (TS = 7) for setting a 95% 
confident-level upper limit on the Fermi-FAT data points. The 
physical interpretation of the SED will be discussed in the next 
section. 


Table 2. Spectral data points of GeV-scale y-ray emission asso¬ 
ciated with RX J1713.7-3946. 


T7 a 

-^mean 

EE{E) 

S[EF{E)} 

7.59 

1^ 

2.40 

0.59 

1.55 

10^ 

4.77 

1.58 

3.14 

10^ 

4.68 

1.71 

6.40 

10^ 

5.94 

1.66 

1.30 

lO'* 

9.08 

2.12 

2.65 

lO'* 

19.0 

3.23 

5.40 

lO'* 

18.6 

4.06 

1.10 

10^ 

9.67 

4.25 

2.24 

10^ 

24.4 

9.60 


Energies are in MeV and the E F{E) flux from RX J1713.7-3946 is 
calculated in units of 10“^ MeV cm“^ s“k 


3. Interpretation 


The har d spectrum of the GeV-band emission initially discov¬ 
ered by Abdo et al. ( [2011| ) and confirmed here suggests an in¬ 
verse Compton (IC) origin of the emission. In a simple one-zone 
scenario for the particle distribution and subsequent emission it 
is difficult to accommodate a strong hadronic component. In this 
section we point out that more complex scenarios are possible, 
although they require some degree of fine-tuning. In an analytic 
assessment and detailed calculations we explore the viability of a 
hadronic-emission scenario involving radiation from the vicinity 
ofRXJ1713.7-3946. 

One-zone models of th e broadband spectra l energy distri¬ 
bution (SED) from X-ray ( [Tanaka et al.||2008') ) to GeV y-ray 
(Fermi-FAT), and TeV y-ray (H.E.S.S.) energies invariably wash 
over, e.g., the variation s in the X-ray spectrum across the rem¬ 
nant ( [Sano et ar|[2014 ), wh ich may reflect electron ageing and 
a changing magnetic field (Rettig & Pohl [2012} Ressler et al. 
[2014[ ). A few general insights can be gained, though. Stringent 
constraints on the gas density in RX J1713.7-3946 are given 
by the absence of a clear detection of thermal X-ray emis¬ 
sion ( [Tanaka et al.[[2008] ), indicating a post-shock gas density 
of ^H^0-8cm“'^. 

Given the low gas density and the y-ray brightness of RX 
J1713.7-3946, hadronic models require a very large energy den¬ 
sity in cosmic rays and a very hard particle spectrum starting 
from lowest energies to reproduce the observed GeV-TeV y-ray 
spectrum. This is difficult to explain by diffusive shock acceler¬ 
ation even in non-linear mode. Feptonic models, on the other 
hand, invariably require an unusually weak magnetic field of 
amplitude ~ 10 yuG, implying that electron energy losses are 
negligible and, at odds with theoretical expectations for efficient 
diffusive shock acceleration, that there is no significant ampli¬ 
fication of magnetic field at the shock. Also, the electron spec¬ 
trum must be moderately soft {a > 2) to simultaneously account 
for X-ray, Fermi and H.E.S.S. data. There are more complicated 
hadronic models, however, which may fit the GeV-scale emis¬ 
sion from RX J1713.7-3946. We shall discuss them below. 


3.1. Hadronic emission from the SNR interior 

RX J1713.7-3946 was likely produced by a core-collapse super¬ 
nova. The progenitor star would have emitted a wind that pushed 
the ambient medium into a dense gas shell at some distance. The 
gas density inside the wind zone is low, implying a weak decel¬ 
eration of the forward shock as is observed. The high velocity 
currently measured suggests that the forward shock has not yet 
reached the wall enclosing the wind-blown cavity. 
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Inoue et aT] ( |2012[ ) suggested that dense gas clouds survive 


passage through the forward shock of RX J1713.7-3946. Slow 
and energy-dependent diffusion of cosmic rays into these clouds 
would then lead t o a hard particle spec t rum a t the high-density 
core of the cloud. [Gabici & Aharoni^ ( |2014| ) demonstrate that 
a fit to the GeV-to-TeV spectrum can be achieved, albeit appar¬ 
ently requiring a mass of ~ 500 Mq inside the SNR, carried 
only by massive clouds that GeV-TeV-scale cosmic rays cannot 
fully penetrate. This is much more than can condense out of the 
wind of the massive progenitor star of RX J1713.7-3946, and 
so the scenario would require that pre-existing clouds have been 
unaffected by the wind. Besides the fact that non-penetration of 
clouds is n ot observed in the y-ray spectra of nearby molecular 
complexes ( [Ackermann et al.|2012 1, we see two difficulties with 
this scenario. 


Firstly, the medium around massive stars tends to be largely 
homogenized before they explode as core-collap se Supernova, 
owing to photo-evaporation and the rocket effect ( [McKee et al.| 
|1984| ). Secondly, as the stellar wind streams around dense gas 
clouds, Kelvin-Helmholtz instabilities followed by Rayleigh- 
Taylor instabilities will disrupt and wash off their outer lay¬ 
ers, leading to large streams of gas with moderate density {n ^ 
1-10 cm“^) that upon contact with the SNR forward shock 
would emit intense X-ray emission. One would have to demon¬ 
strate that its intensity is compatible with the published upper 
limits for thermal X radiation. 


A large fraction of the gas clouds and the outer gas shell 
will thus be located outside of the SNR, where they would be 
illuminated by the cosmic-ray precursor to the forward shock. 
This scenario we shall discuss in the next subsection. 


3.2. Hadronic emission from the cosmic-ray precursor 

The width of the cosmic-ray precursor to astrophysical shocks 
is strongly energy dependent on account of the momentum de¬ 
pendence of the mean free path for scattering. At a given loca¬ 
tion, this leads to very hard particle spectra with cut off at an 
energy determined by the distance from the shock, the diffusion 
coefficient, and the shock velocity (Gabici et al.|2009t[Ellison &| 
[Bykov[2QTT|[Telezhinsky et al.[2012b[ ). It has been argued before 

that the escape of cosmic rays from SNRs might lead to strong 
TeV-band y-ray emission from molecular clouds located 10 pc 
or more from the remnant. In the case of RX J1713.7-3946, the 
target gas has to be located very close to the forward shock, 
otherwise we cannot observe intense emission at 50-500 GeV. 
Consequently, it is illuminated not primarily by so-called run¬ 
away cosmic rays that have escaped from the system, but by 
the cosmic-ray precursor to the forward shock. We do not doubt 
that it would be possible to find parameters that permit repro¬ 
ducing the observed y-ray spectrum of RX J1713.7-3946 within 
such a scenario. Rather, we predict a gradual but detectable in¬ 
crease in y-ray flux over a decade or two that renders this sce¬ 
nario testable. To substantiate our claim, we shall first present 
an analytic estimate to illustrate why the signal would grow with 
time. We will then use a more sophisticated model to calculate 
the time dependence of the y-ray spectrum. 

3.2.1. Analytic estimate 

We shall begin with an analytic model that despite its simplic¬ 
ity will demonstrate the existence and magnitude of time de¬ 
pendence in the y-ray signal from a gas cloud just outside of 
an SNR. The transport equation for the differential number den¬ 


sity of cosmic-rays, N(E, r), comprises terms for spatial diffu¬ 
sion and advection, that assuming spherical symmetry can be 
written as 


dN(E, r) 1 d 


dt 


dr 


2 n,p .dN{E,r) 
r D{E, r) 


dr 


■r^ V(r)N(E, r) 


^DEN(E,r) = Q(E,r), 


(4) 


where T>e denotes the differential operator in energy. Rewriting 
Eq.j^in co-moving coordinates anchored on the shock position, 
rs, 

r ^ p = r-r^ and Vp = V - r^ (5) 

and restricting to the vicinity of the forward shock (p r^), 
we may assume a plane-parallel geometry. As particle accelera¬ 
tion is faster than the hydrodynamical evolution of the SNR, the 
particle distribution in co-moving coordinates can be in good ap¬ 
proximation written in the steady-state limit. 


d 

dp 


D(E,p) 


dN(E,p) 

dp 


VpN(E,p) 


^DeN(E,p) = Q(E,p). 


( 6 ) 

Note that in the shock precursor region sources are absent, i.e. 
Q{E,p) = 0. The term containing De is expected to be small, 
because energy-changing processes operating in the precursor 
region are very slow compared to the SNR evolution, and so we 
have to search for homogeneous solutions to the spatial part of 
the transport equation. 


dp 


dN(E,p) 

D(E,p) ; -VpN(E,p) 
dp 


0 . 


(7) 


Generally, there are 2 solutions to Eq.j^that are approximately 


^ N 2 -C exp 

r 0 




Vn 


D(E,p')} ’ 


( 8 ) 


where A and C are constants that are determined by matching 
the cosmic-ray density at the shock. The plasma flow speed in 
the precursor region can be somewhat modified by non-linear 
shock modification, and so we only know Vp > -/*§. It is obvi¬ 
ous that only N 2 satisfies the boundary condition N{p = 00 ) = 0, 
therefore it must describe the true cosmic-ray density profile rea¬ 
sonably well. 

Now suppose a gas cloud is located at radius ro, corre¬ 
sponding to the time-dependent co-moving coordinate po(0 = 
^0 “ ^s(0 - ^0 “ where the last expression is approximately 
valid for short time intervals. The gas cloud is illuminated with 
cosmic rays of differential density 


( 9 ) 

The time scale on which the cosmic-ray illumination changes 
can be estimated as 

_ N _ D(E,po) _ D(E,po) 

^ If I IPoVpl • 

If we scale the diffusion coefficient to the Bohm limit, D = 
7] r^ c/3, wc find for the time scale of variation 


(10 yrs) 77 


(- 


E \ 


B 


110 TeV/ \ lOpG/ \ 3000km/s 


( 11 ) 
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Protons with 10 TeV in kinetic energy produce y-rays at energies 
around 500 GeV, at which imaging atmospheric Cherenkov tele¬ 
scope have their peak sensitivity. As long as the diffusion coef¬ 
ficient is not significantly larger than the Bohm limit, we expect 
observable variations in the y-ray flux from bright, young SNRs 
such as RX J1713.7-3946, in particular if streaming instabilities 
amplify the magnetic field to amplitudes in excess of 10 yuG. 

To be noted is that the variation time scale does not explic¬ 
itly depend on the distance from the shock. Implicitly, a certain 
distance dependence arises from the fact that the diffusion coef¬ 
ficient must be evaluated near the location of the gas cloud, not 
at the shock. 


3.2.2. Detailed modeling 


In this subsection we describe a more detailed model of RX 
J1713.7-3946 and its high-energy (HE) y-ray flux variability 
that we expect according to the predictions outlined in subsec¬ 
tion |3.2.1[ We consider a system consisting of a core-collapse 
SNR and a shell of dense gas located just outside the remnant. 
We do not pretend to give a good broadband fit to the data be¬ 
cause that would require a sophisticated model of the evolu¬ 
tion of, and subsequent emission from, the remnant itself (e.g. 
[Telezhinsky et al.||2013| ). In any case, the remnant itself is a 
rather faint source of y-ray emission on account of the low- 
density environment left by the progenitor stellar wind. The 
hadronic emission is dominated by y-rays from the dense shell 
that is illuminated by high-energy cosmic-rays in the precursor. 
Consequently, the overall spectrum is expected to be hard. In 
low-energy (LE) y-rays the flux from the SNR and the shell may 
be comparable, because a low gas density near the forward shock 
is compensated by a high intensity of LE cosmic rays. The flux 
of HE y-rays from the dense shell and its variation depend on 
the cosmic-ray diffusion coefficient and the current kinematics 
of the shock, parameters that change little if sophisticated mod¬ 
els of SNR are considered. Therefore, the prediction of HE y-ray 
flux variation offered here is robust. 

We numerically solve the full transport equation for cosmic 
rays 0 on a grid that extends from the SNR center to several 
dozens of SNR radii ahead of the forward shock (ES). The grid 
extension allows to account for the escape from the precursor 
of particles of both species, protons and electrons. We can use 
either full hydrodynamical simulations or analytic solutions to 
trace the evolution of the SNR, i.e. the velocity and radius of 
shocks and the plasma flow profile inside the remnant. A crucial 
parameter for particle acceleration is the diffusion coefficient. 
We assume it to be Bohm-like, Db = rjrL v/3, where rj = 1.0 (at 
the shock), ri is the particle Larmor radius, and v is the parti¬ 
cle velocity. Ahead of the shock the diffusion coefficient linearly 
increases ( Bell| 197 8} [Drury || 1 983| ) until thirty SNR radii ahead 
of the shock it reaches Dq = 10^ (E/10 GeV)^/^ , a 

value typical for the galactic propagation. Close to the shock dif¬ 
fusion is still Bohm-like, and so a typical cosmic-ray precursor 
will be established. Eurther ahead, the probability of escape from 
the system becomes significant, and particles may leave the sys¬ 
tem. Eor full details we refer a reader to a series of publications 
( [Telezhinsky et al.|2012a|b[[2013| ). 

Given that thermal X-rays are not detected, the current num¬ 
ber density of gas upstream of the blast wave must be around 
0.05 cm“^, much lower than the typical density in the interstellar 
medium. Therefore, the explosion should have taken place in the 
wind-blown bubble of the massive progenitor star, and the SN 
may be classified as a core-collapse event. To reconstruct these 


parameters we used analytical solutions for the free-expansion 


stage of SNR expanding into a wind-blown cavity (Truelove & 
McKee||1999| ). The current SNR radius R ^ 8.65 pc and shock 


speed Vsh - 4200 km/s (for a distance of 1 kpc) permit a range 
of parameter values for the mass-loss rate and the progenitor 
wind speed. As example we here assumed for the mass-loss rate 
M ^ 10“"^ Mq, for the wind speed ^ 50 km/s, and for the 
explosion energy Esn = 10^^ ergs, corresponding to a relatively 
fast red-supergiant wind. It serves as demonstration of the sta¬ 
bility of our results against the choice of hydrodynamical pa¬ 
rameters that we even tested a Sedov solution, and the implica¬ 
tions for the fiux increase at a few 100 GeV were virtually the 
same. The magnetic field in the precursor region is assumed to 
be Bp = 23yuG which is shock-compressed to a downstream field 
Bd = 75yuG inside the SNR ( |Ballet|2006}|Acero et al.|2009| ). 

RX J1713.7-3946 is considered to be 1600 years old, and 
we use that as fiducial age of RX J1713.7-3946. Eitting the 
Truelove-McKee profiles permits a certain liberty in choosing 
the distance, the age, and the current shock velocity. Note that 
the actual age is of minor relevance for the scenario at hand, 
important are the time span between 2 observations of the ob¬ 
ject and, as Eq. illustrates, the shock velocity and the diffu¬ 
sion coefficient. The allowed range in current shock velocity is 
fairly small, thus setting limits on the age, and despite entering 
Eq. p^quadratically it does not lead to a large uncertainty in the 
timescale of fiux variability. The location of the dense gas shell is 
prescribed by the necessity to reproduce the presently observed 
y-ray spectrum from RX J1713.7-3946, as the precursor of cos¬ 
mic rays at 10 TeV should just reach it. That leaves the diffusion 
coefficient and in particular the efficacy of magnetic-field am¬ 
plification as the main sources of uncertainty, for which we use 
conservative assumptions. 

We numerically solve the time-dependent cosmic-ray trans- 
port equation (Eq. using the above parameters ( jTelezhins^ 


et al.|20T2^ . We then calculate emission spectra and inten¬ 


sity maps at two epochs, at 1600 and at 1610 years, including 
pion-decay ( [Huang et al.|2007| ), inverse Com pton (IC) and syn - 
chrotron emission in turbulent magnetic field ( Pohl et al.|2014[ ). 
The dense gas shell radially extends from 8.75 pc to 9.05 pc from 
the SNR center, and so the shock does not interact with the dense 
material, but CRs in the shock precursor can efficiently reach the 
dense gas. A small separation between the gas shell and the for¬ 
ward shock is necessary to reproduce the v Ey-peak at ~ 1 TeV 
in the hadronic y-ray emission. The results of our calculations 
are presented at Eig. (solid red line). The HE y-ray spectrum 
at the age of 1600 years is normalized to fit H.E.S.S. data and 
to reconstruct the measured flux, F = (17.2 + 0.4) x 10“^^ pho¬ 
tons cm“^ s“^ ( Aharonian et al.|2()TT I. 

The parameters are defined by the gamma-ray fiux ratio be¬ 
tween the TeV band and the GeV band as well as the location 
of the peak in the vFy spectrum. The latter determines the sep¬ 
aration between the forward shock and the gas shell, somewhat 
depending on the scaling of the diffusion coefficient in the pre¬ 
cursor. The former prescribes the mass or gas density in the shell, 
~ 300 Mq and nu ^ 35 cm“^ for our SNR expansion model with 
17 Mq of swept-up gas in the remnant. This is simply the amount 
of target material needed to fit the observed gamma-ray spec¬ 
trum; a thorough study of the hydrodynamical circumstances 
that could give rise to a shell with the required mass is beyond 
the scope of this paper. The IC contribution is constrained by the 
observed fiux of synchrotron X-rays and our choice of magnetic- 
field strength. Reproducing the spectral width of the y-ray signal 
in the vFy spectrum shown in Eig. [^requires a diffusion coeffi¬ 
cient well in excess of the Bohm limit. In other words, the pa- 
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Fig. 6. Broadband y-ray spectrum of RX J1713.7-3946 as mea¬ 
sured with LAT (our results, black points) and with H.E.S.S. 
(blue points), compared to calculated HE y-ray spectra at the age 
of 1600 years (solid red line) and 1610 years (dashed red line). 
The IC contribution to both spectra is indicated for additional 
information (gray dotted line). The relative difference between 
spectra at two epochs is shown in the lower panel (solid green 
line). 

rameter rj ^ 10 in Eq. [TT] increasing with distance from the 
shock, and so with B = 23 jjG and = 4200 km/s we expect as 
time scale of variation r ^ 100 years. The numerical treatment 
indicates a flux increase of 15% over ten years, consistent with 
the analytical estimate. 

The main prediction of our model thus is a secular increase 
with time of the HE y-ray flux that is largely independent of the 
choice of parameters, provided they permit reproducing the peak 
in the vFy spectrum at ~ 1 TeV for the current size and expan¬ 
sion rate of RX J1713.7-3946. Eig. (dashed red line) shows 
the evolution of the y-ray over 10 years. The relative increase 
is about 15% at around 500 GeV (see the green line at Eig. [^. 
Note that the IC flux will not change during this period because 
the amount of target photons interacting with electrons in the 
shock precursor does not change noticeably. Also note that the 
change in morphology will be too small to be detectable with 
CTA. 

The calculated flux in the 1-100 TeV energy band at the age 
of 1600 years amounts to Fi6oo = 17.2x 10“^^ photons cm“^ s“\ 
compatible with that measured with H.E.S.S., whereas 10 years 
later the flux would beFi6io = 19.9x10“^^ photons cm“^ s“\ 
with most of the increase coming from the shell. This change 
in the flux could be potentially detected with the CTA ( |Acharya| 
|et al.|2013] ). The systematic uncertainties in the flux scale for the 
original measurement with H.E.S.S. are probably too large for 
that, and so one would need CTA data taken at least 10 years 
apart. 

4. Conclusions 

In this work we analyzed more than 5 years of Eermi-LAT ob¬ 
servations of the SNR RX J1713.7-3946, using high-resolution 
y-ray background maps that account for the HI self-absorption. 
RX J1713.7-3946 appears as circular source with radius r ^ 
0.55° that is detected at 13cr level in the energy range between 
500 MeV and 300 GeV. We see a tendency in the GeV band 


for the same intensity enhancement in the northwest region of 
the remnant that is seen with H.E.S.S. and probably due to par¬ 
ticles interacting with denser gas located there. The measured 
spectrum of RX J1713.7-3946 indeed shows a hard spectral in¬ 
dex ofy= 1.53 + 0.07 with an integral flux above 500 MeV of 
F = (5.5 ± 1.1) X 10“^ photons cm“^ s“^ 

Phenomenological fits of the broadband SED demonstrate 
that in simple one-zone models for particle acceleration and 
emission it is hard to accomodate any pion-decay scenario. 
These fits also show that a common particle power-law spectrum 
with an exponential cutoff for electrons and protons faces some 
difficulty in simultaneously accounting for the X-ray, GeV y-ray, 
and TeV y-ray observations. A more complex model with addi¬ 
tional parameters and more flexible assumptions could modify 
the shape of the broadband spectrum and give a better descrip¬ 
tion of the y-ray data. 

We discussed alternative hadronic scenarios and pointed out 
some difficulties with them. In particular, we introduced a model 
of hadronic emission from RX J1713.7-3946 that reproduces the 
hard GeV spectrum with index y ^ 1.5 as emission from a shell 
of dense gas that is located a short distance upstream of the for¬ 
ward shock. The existence of such a shell may be related to the 
winds of the progenitor star. Although this hadronic scenario 
is somewhat speculative, it has the benefit of being testable: 
We predict the hadronic y-ray flux to increase at the level of 
15% over 10 years which should be observable with the future 
CTA facility that is p rojected to operate for more than 10 years 
( Acharya et al.|2013] ). We have verified that the growth in y-ray 
flux is largely independent of the choice of parameters, provided 
they reproduce a v Fy-peak at ~ 1 TeV Whereas our model as¬ 
sumes spherical symmetry, real SNR are not that simple, and 
conditions are not the same in different regions. The stability of 
our result suggests that such a variation of parameters does not 
significantly reduce the secular increase in y-ray flux. 

A second implication of the scenario is that, within the next 
50 years, the forward shock of RX J1713.7-3946 should hit the 
gas shell, which would then light up in X rays and permit un¬ 
precedented studies of shock-cloud interactions. We stress that 
similar intensity fluctuations are expected in models invoking 
the forward shock moving through a medium in which most of 
the mass is organized in dense gas clouds. 
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